Spontaneous formation of a tt soliton in a superconducting wire with an odd number 
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We consider a one-dimensional superconducting wire where the total number of electrons can be 
controlled in the Coulomb blockade regime. We predict that a tt soliton (kink) will spontaneously 
form in the system when the number of electron is odd, because this configuration has a lower 
energy. If the wire with an odd number of electrons is closed in a ring, the phase difference on 
the two sides of the soliton will generate a supercurrent detectable by SQUID. The two degenerate 
states with the current flowing clockwise or counterclockwise can be utilized as a qubit. 

PACS numbers: 74.50.+r 03.67.-a 05.45.Yv 73.23.Ra 
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Introduction. In a mesoscopic sample in the Coulomb 
blockage regime, the total number of electrons can be 
controlled discretely by gate voltage . It was predicted 
theoretically |2j and verified experimentally |^ that the 
ground-state energies of a superconducting sample with 
even and odd number of electrons differ by the supercon- 
ducting gap Aq, because the odd electron cannot form 
a Cooper pair. In the present paper, we discuss an- 
other manifestation of the even-odd asymmetry in a one- 
dimensional (ID) superconducting wire or constriction. 
We show that when such a system has an odd number of 
electrons, it can lower its energy by spontaneously form- 
ing a kink soliton of the order parameter, analogously to 
the charge-density-wave solitons in polyacetylene 
and other systems with similar broken symmetry 
the wire is closed in a ring, the phase difference on the two 
sides of the soliton would produce a supercurrent in the 
ring, which could be detected by a sensitive SQUID. The 
two degenerate states with the current flowing clockwise 
or counterclockwise can be considered as a qubit . 

ID superconductor with a kink soliton. Let us con- 
sider a ID s-wave BCS superconductor characterized by 
the mean-field Hamiltonian density 




H{x) c| (x) i C| (x) + c| (x) i ci (x) 
+ ci (x) A* {x) ct {x) + c| {x) A{x) c| (x) 



(1) 



|A(x)|V5. 



Here cj.^|(a;) are the creation and destruction operators 
of electrons with the spin projections t and J, at the point 
X, A(x) is the superconducting pairing potential satisfy- 
ing the self-consistency condition A(a;) = g{ci{x)ci{x)), 
5 < is the BCS coupling constant, and ^ is the kinetic 
energy operator. Near the Fermi energy, the electron dis- 
persion relation can be linearized: ^ « —aivpdx, where 
a = ± labels the right- and left-moving electrons with 
the momenta near the two Fermi points akp, vp '^s the 
Fermi velocity, and the Planck constant is ?i = 1. All 
calculations in the paper are done at zero temperature. 

The eigenstates of Hamiltonian (|l|) are the Bogoli- 
ubov quasiparticles characterized by the wave functions 



'ijj^n\x) — [u\'i'\x),v^n\x)] that satisfy the Bogoliubov- 
de Gennes (BdG) equation 



(2) 



with the eigenenergies The spectrum of Eq. (|^) 

consists of solutions with positive and negative ener- 
gies related by charge conjugation. The superconducting 
ground-state energy F relative to the normal state energy 
and the BCS self-consistency relation can be expressed 
in terms of these solutions [01: 



F 



Sri, a (-^" ^ + IC" ^ 

A(a;) = -gE 



/dx|A(x)|V5, (3) 
(4) 



(a) ^ N (a)*/ \ 
Uii [X) Vn [X), 



where the sums are taken over the occupied states with 
< 0. 

Suppose the pairing potential has a phase soliton 
(kink) at X = 0, so that A(x) Aqc^'-'^ a.t x ±oo: 

A{x) = Ai+iA2{x), Ai = Aocos0, (5) 
A2(a;) = KWi? tanhKx, k = Aq sin(/)/wi?. (6) 

In this case |ll[ , the spectrum of Eq. (||) has a contin- 
uous part ±Eq°'\ where q is the quasiparticle momentum 
counted from akp, and discrete subgap states Eq"^ : 

= ^J{vpqr + AI = -aAo cos 0. (7) 

The energies E^{(j>) of the subgap states (|^) are shown 
in Fig. |l|(a) by the thick solid and dashed lines. The to- 
tal energy F, given by Eq. (||), can be separated into the 
contributions from the continuous spectrum and from the 
discrete states. The former contribution, counted from 
the ground state energy of a uniform BCS superconduc- 
tor with an even number of electrons, has the form 



F, = Ao-2j2{E,^-E,^)- J 



dx 



\A{xr-Al 



(8) 
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FIG. 1: (a) The energies -Eg"' of the subgap bound states vs 
phase (/!>. The thick solid and dashed hnes correspond to the 
transmission coefficient t — 1, and the thin lines to r = 0.9. 
(b) The total energy F of the one-channel superconducting 
wire with an odd (solid line) and an even (dashed line) number 
of electrons vs phase (j>, calculated for r = 1. 



The first term in the r.h.s. of Eq. (J^ reflects the fact that 
one state has been removed from the continuum to the 
discrete spectrum. The factor 2 in front of the second 
term represents the equal contributions from the states 
with a = ±. In the second term, qj — qj — Cj / L is the 
quasiparticle momentum displaced by the phase shift Q 
from the unperturbed value qj = 27rj/L for a uniform 
system without a kink. (Here the periodic boundary con- 
dition over the length L is implied, and j is an integer.) 
The calculations of the phase shifts and the dependence 
of Fc on ip are given at the end of the paper (see also 
|, |, H). The result is 



2Ao 



(9) 



Let us examine how the ground state energies Fe((j)) 
and Fo{4>) of the system with an even and an odd num- 
ber of electrons depend on the phase difference 2(f). At 
= 0, the spectrum has only the continuous part, which 
is completely filled below — Aq and completely empty 
above -I-Aq when the number of electrons is even. As 
(/) increases, the discrete state (|^) with the negative en- 
ergy remains occupied, whereas the one with the positive 
energy is empty. Thus, 



Fe((/)) =i^c(<^)- Ao|cOS( 



(10) 



As shown in Fig. |l|(b) by the dashed line, Fe{(t>) is mini- 
mal at = 0. Thus, the system with an even number of 
electrons prefers a uniform state. 

On the other hand, when the number of electrons is 
odd, the quasiparticle state with the energy -I-Aq is filled 
by the odd electron at = 0. Thus, as increases, 



both the lower and upper discrete states (^ stay occu- 
pied. Because their energies cancel, the total energy in 
the odd case is 



(11) 



which is shown by the solid line in Fig. |^(b). In con- 
trast to the even case, Fo{<t>) decreases with increasing 
(j) and has a minimum at = 7r/2 with the energy 
Es — Fo{tt/2) — 2Ao/7r < Aq, lower than in the uniform 
state. Thus, a ID superconductor with an odd number 
of electrons should spontaneously create a kink with the 
phase difference 20 = tt in order to decrease the total 
energy. This phenomenon is analogous to creation of a 
soliton in a charge- or spin-density wave with an odd 
number of electrons H, ^, 0| . 

Now let us consider the case where the ID supercon- 
ducting wire has N channels, e.g. consists of N parallel 
chains with the common potential A(a;). When the num- 
ber of electrons is odd, N discrete states with negative 
energies (0) and one with positive energy are filled. So, 
the total energy of the system is Fo{4>) = NFc{(j)) — {N — 
l)Ao|cos0|, and the minimum (Ao2iV/7r) sin(7r/2iV) is 
achieved at 20 = tt/N. Thus, the soliton also forms in 
this case, but with the phase difference reduced by 1/iV. 
When TV ^ 1, the energy gain is insignificant, and soliton 
formation becomes irrelevant. 

Superconducting ring with a kink. Now suppose that 
the ID superconducting wire (with iV = 1) is closed in 
a ring of length L. The length L is assumed to be much 
longer than the coherence length ^ — hvp /Aq, so the 
soliton energy Eg is not significantly affected by closing 
of the wire. On the other hand, L should be short enough, 
so that the whole ring is in the Coulomb blockade regime, 
and an external gate can control the total number of 
electrons to be odd. The latter condition requires that 
Ao < i?c 111, where Ec ^ e'^/D is the Coulomb charging 
energy, and D is the overall size of the sample. For a 
simple ring, D = L, but in general Z? > L, as discussed 
below. Thus, we need 



i = Tivf/Ao <^L<D< eVAo, 



(12) 



so the Fermi velocity has to be much lower than the Bohr 
velocity: vp ^ /% — 2.2 x 10* cm/s. This condition is 
satisfied in a typical metal with vp ~ 10^ cm/s. 

When the number of electrons is odd, the soliton cre- 
ates the phase difference 20 = tt between the ends of 
an open wire. However, when the wire is closed, the 
superconducting phase 9 must continuously interpolate 
between the phases ±7r/2 at the opposite sides of the 
soliton, so A(a;) = Aoe'^*^^-'. The small phase gradient 
\d9/dx\ — -k/L creates the supercurrent 



J 



^erig do 
2m dx 



± 



evp 



(13) 



where rig = Akp/2'!T is the superfluid density equal to the 
total electron density at zero temperature p^. Notice 
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that the supercurrent J can flow in the ring either clock- 
wise or counterclockwise, thus the time-reversal symme- 
try is spontaneously broken. These two degenerate states 
can be viewed as a qubit The supercurrent ( p^ is 
equal to the maximal persistent current in a normal (non- 
superconducting) ring with an odd number of electrons 
iQ . It corresponds to the odd electron of charge e mov- 
ing with the Fermi velocity vp around the ring of cir- 
cumference L. This current creates the magnetic flux 
$ ^ JL/c ~ evp/c ~ 10~^$o (where $o — Trhc/e is the 
flux quantum), which could be detected, in principle, by 
a sensitive SQUID. We predict that, when the number of 
electrons in the ring is switched between odd and even by 
an external gate, the current ( |7^ should appear and dis- 
appear in the ring. In the case of A'^ channels, the current 
is given by the same Eq. ([l3|). 

Above we considered an idealized situation of a uni- 
form superconducting ring. However, the conclusion also 
applies to other, more realistic experimental setups. Sup- 
pose a straight ID superconducting wire of length L, e.g. 
a carbon nanotube [p^ , is shunted by a thick conven- 
tional superconductor of size D. In this case, the soliton 
and the phase gradient are localized within the ID wire, 
whereas the phase gradient over the shunt is close to zero, 
so the supercurrent is given by Eq. (p^. Another setup 
is a narrow constriction with a small number of con- 
ducting channels in a thick superconducting ring. In this 
case, L is the length of the ID constriction, and D is 
the overall size of the ring. In general, when the ring 
is nonuniform, the transmission coefficient t through the 
weak link is less than one. In this case, the energies of 
the subgap states are Eq = ±Ao\/l — rsin^ pT| , 
shown by the thin lines in Fig. |^(a). Although we cannot 
calculate Fc{(j)) exactly when r 7^ 1, the total energy of 
the system should be close to Eq. (^ if r « 1, and the 
kink soliton should be still energetically favorable. Thus, 
our theory should be applicable to the atomic contacts 
with N ^ 1 and r « 1 realized in experiments p^ ]. 

Calculation of Fc{(j)) via phase shifts. To simplify pre- 
sentation, we write the equations below specifically for 
a = +. In order to calculate the phase shifts, let us in- 
troduce the wave functions f±{x) — u{x) ± v{x). (Notice 
that this index ± is different from a.) Then the BdG 
equation (^ acquires the supersymmetric form |p7| , 



[E T Ai)/± = ^ivF df^/dx T ^^2{x)f^, (14) 
dA2{x) 
dx^ dx 



E'+vl-i^-\Aix)\'Tvp 



/±-0. (15) 



Substituting Eqs. (|) and (|) into Eqs. (|T|) and (|l|), we 
find the (unnormalized) wave functions 



/+(x)=e*''^ /_(x) 



vpq + iA2ix) ^^^^ 
E + Ai 



From Eq. (|T|), we observe that, when x changes from 
—00 to +00, /_ acquires the phase factor S = (g -|- 



ik) / {q — in). In order to impose periodic boundary condi- 
tions consistently, let us assume that the system also has 
an antikink far away from x = 0. In a similar manner, it 
is easy to show that /+ acquires the same phase factor 
S upon traversing the antikink. Thus, when x changes 
from one boundary of the sample to another, both f± ex- 
perience the phase shift C,q = arctan(K/g) per one soliton 
|l^ . The same result can be also obtained by imposing 
zero boundary conditions at the ends of the wire. 

Using the relation qj — qj — Q/L, we write Eq. (^) as 



Fc = Ao 



2KVp 

9 



(17) 



where we took the integral over x using Eqs. (g) and (|g). 
Integrating by parts over q and taking into account the 
BCS self-consistency relation 



dq 

2'KEn 



(18) 



we rewrite Eq. (17) as follows 



F, 



Substituting Qq — arctan(K/g) into Eq. ( p^ and taking 
the integral over g, as in Ref. we obtain Eq. (||). 

An alternative calculation of Fc{(j)). Following the 
method of Ref. B , let us calculate the derivative 



dF,(0) 



dx 



SFc dA*{x) 



(20) 



_SA*{x) 84) 

We write the eigenenergies of the BdG equation (|^) as 

= j dx [4:^>{x)iui^\x)-v^^\x)iv\^>{x) 
+ zil") {x) A* (x) vi"> (x) + (x) A{x) ui"> (x)], (21) 



where the normalization condition J dx{\u\'^ + \'i 
is implied. Taking variational derivatives of Eqs. 
(pT|) with respect to A* (a;), we find 



and 



SFr, 



SA*ix) 



5:4")(x).(")*(x)-^. (22) 



q,a 



The normalized wave functions ^pq 
deduced from Eq. ( p^ ) 



q , Uq 



' Eq + aAi 
^LE„ 



£,+aAi 



) can be 



(23) 



(16) Using Eq. (H), we find 



5^<)«(")* = A(x) j 

q,a 



dq 



k{t:/2 ~ (j)) 



2TTEq 27rcosh^(Ka;) 



(24) 
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where the summation over qj has been simply replaced 
by integration over g, because the phase shifts contribute 
only a negligible term proportional to 1/L. Substituting 
Eq. (p|) into Eq. (|2|) and using Eq. (Ill), wc find 



(25) 



(5A*(x) ~ 27rcosh^( KX) 



Substituting Eq. ( p5| ) and OA* /d(j) determined from Eqs. 
(^ and (^) into Eq. (^0|) and integrating over x, we get 



dFr 



2Ao /TT A . , 



(26) 



Integrating Eq. (26) over (f>, we recover Eq. (^. 

We can also check that the self-consistency condition 
(^ is satisfied a.t — tt/2. Indeed, the last term in 
the r.h.s. of Eq. (^J) vanishes, and the first term gives 
—A{x)/g because of Eq. (|l8|), whereas the contributions 
of the localized states 



^0 ^ ' " 2coshAtx I -a 



(27) 



mutually cancel for a — ±. 

Conclusions. We have shown that a ID supercon- 
ducting wire with an odd number of electrons can lower 
its energy by creating a tt soliton (kink) of the order 
parameter, so that the odd unpaired electron occupies 
the midgap bound state localized at the soliton. In or- 
der to exhibit this effect, the wire does not have to be 
completely uniform. It is sufficient that it has a ID por- 
tion with the length L longer than the coherence length 
^ and the atomic-size cross-section with the number of 
conducting channels A'^ ~ 1 and the transmission coeffi- 
cient r « 1. Such atomic Josephson contacts have been 
already realized experimentally Eq] . Another possibility 



is to use a carbon nanotube 1 15 



If the wire is closed in a ring, the phase difference on 
the two sides of the soliton would generate a supercur- 
rent in the ring, which could be detected by a sensitive 
SQUID. The total number of electrons in the ring can be 
controlled by an external gate in the Coulomb blockade 
regime. We predict that the supercurrent should appear 
in the ring when the number of electrons is odd and dis- 
appear when the number is even. The current can flow 
either clockwise or counterclockwise. These two degen- 
erate states can be viewed as a qubit and utilized for 
quantum computing 

Spontaneous current in a ring with phase circulation tt 
was recently predicted for unconventional superconduc- 
tors using exotic theory of fractional spin-charge sepa- 
ration Q and was disproved experimentally |Q. That 
theory did not consider the difference between even and 
odd number of electrons or the setup with only one or 
few conducting channels. These mesoscopic aspects play 
crucial role in our theory, which is formulated completely 
within the conventional BCS framework. 
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